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Abstract 

We study the behavior of the fermion determinant at finite tem- 
perature and chemical potential, as a function of the Polyakov loop. 
The phase of the determinant is correlated with the imaginary part 
of the Polyakov loop. This correlation and its consequences are con- 
sidered in static QCD, in a toy model of free quarks in a constant A 
background, and in simulations constraining the imaginary part of the 
Polyakov loop to zero. 

1 Introduction 

The experimental study of QCD at finite density and temperature, and the 
possible observation of the quark-gluon plasma are the focus of activity of an 
increasing fraction of the physics community. A parallel study by computer 
simulations would be highly desirable. However, after integration over the 
fermions, the integrand p(U) in the finite density partition function is com- 
plex. The measure used in Monte Carlo simulations is (..)mc = d[U] \p(U)\, 
and observables are computed as (O) = (O p{U) /\p(U)\) mc I (p(U) /\p(U)\) mc ■ 
The denominator is exponentially small in the 4-volume of the lattice, and 
measuring it with a given accuracy requires an exponentially growing amount 
of statistics. This is the so-called sign problem. Over almost twenty years, 
progress in controlling the sign problem has been negligible. 

For this reason one has turned to approximations to QCD: the strong 
coupling limit, and more recently the static limit of QCD have attracted 
a lot of interest jl], [|]. Static QCD is still affected by the sign problem, 
which limits investigations to small volumes, far from the continuum limit. 
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It would be very desirable to extend the domain of study of static QCD to 
more realistic regions of parameters, and to go beyond the static approxi- 
mation and allow for at least some spatial hopping of the quarks. Our paper 
proposes and tests a method to do this. 

We have analyzed static QCD results of |l| and found a strong correlation 
between the phase of the determinant </>(det) and the imaginary part Pi of 
the Polyakov loop P. We study this correlation in Section I, and show that 
it is to be expected for heavy quarks. 

This provides the motivation for a simple toy model, studied in Section 
II, where all degrees of freedom are suppressed except the Polyakov loop: 
free quarks in a gauge field background with constant temporal component 
Aq and vanishing spatial components. This model can be solved easily in 
Fourier space. It shows several features of QCD. In particular, it exhibits the 
ordering effect of the chemical potential /i, similar to that of the temperature 
T. Under simple assumptions, a phase diagram ^ C (T) can be obtained at 
small [i. A mechanism for a transition at low temperature is proposed. 

Finally, the correlation between the phase of the determinant and the 
imaginary part of the Polyakov loop can be put to profit by performing 
constrained simulations, where Pi = 0. The sign problem is reduced. The 
validity of imposing such a constraint, and the benefits this constraint pro- 
vides, are reviewed in Section III. 
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Figure 1: Correlation between <^>(det) and Pi in static QCD, on a lattice 
Ll x 2 (L s = 6, 8, 10, C = 1/20, j3 = 5.0). On the right, Pi has been rescaled 
by \/Lj and 0(det) by 1/i/Zf. The data are from Ref.[l]. 
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Figure 2: Same as Fig. 1 with C = 1/0.6 and j3 = 5.0 on a 6 3 x 4 lattice. 

2 Correlation between Pi and </>(det) 

We have reanalyzed the static QCD data of Recall that static QCD is 
obtained by taking the quark mass m q and the chemical potential simulta- 
neously to infinity, such that the quark density is finite. This is achieved 
by keeping constant the parameter C _1 = (e^ /2m q ) Nt , where N t is the 
temporal extent of the lattice. The quarks are static, but the fermionic 
determinant is still sensitive to changes in the gauge field and the quark 
density. Fig. la shows almost perfect correlation between the phase of the 
determinant and the imaginary part of the Polyakov loop, for 3 spatial vol- 
umes V. Furthermore, as V is increased, one observes that the fluctuations 
in Pi vary like l/y/V, as any extensive variable, and that fluctuations in 
c/>(det) increase like yV . Under this rescaling, the results for the 3 different 
volumes become indistinguishable (Fig. lb). 

As the chemical potential (equivalently C in static QCD) increases, or 
as the temperature decreases, fluctuations in the phase increase, and the 
correlation becomes harder to unravel because of the compact nature of the 
phase. Fig. 2 shows the worst case studied by Blum et al. 0. If one defines 
0(det) as 5Zj<KAj), with 0(Aj) g] — it, +7r], where Aj are the eigenvalues of 
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the Dirac matrix, the correlation remains significant. 

This correlation should be expected. In the loop expansion of the fermion 
determinant, log(det(l — kM) = —S e ff = TrM*, where k ~ l/m q , 

the first loop affected by the chemical potential fi is the (timelike) Polyakov 
loop. The corresponding contribution to S e fj is 

2(2k) n * J2(cosh(N t iJ,)P r (x) + i smh(N t fi)P l (x)) (I) 

X 

i.e. (f)(det) = 2(2K) Nt V sinh(N t n)Pi , (2) 

where the summation J2x ^ s over au spatial coordinates x. Higher-order 
terms spoil this perfect correlation. They become more important as Nt 
increases or m q decreases. In static QCD, there is no spatial hopping of the 
quarks, so that the only additional terms are powers of the Polyakov loop 
at each site, which are suppressed by powers of2(2 K ) 7Vt sinh(iV t /z). 

3 Toy model: free quarks in a constant Aq back- 
ground 

We have just seen that static QCD is in fact a Polyakov loop model, where 
the phase of the determinant is almost equal to the imaginary part of the 
Polyakov loop. But even in static QCD, calculations are still severely lim- 
ited by the sign problem. Monte Carlo integration performs an averaging 
over the gauge fluctuations and the fluctuations of the Polyakov loop P(x) 
at each spatial point x. We can simplify this model further by keeping only 
the zero-mode of such fluctuations and neglecting, for now, the other fluctu- 
ations of the gauge field. We end up with a single SU(3) degree of freedom, 
the Polyakov loop P. Our model consists of free quarks in a gauge field back- 
ground with zero spatial components and constant temporal component Aq. 
Hence, P = exp(iA^aAo), where a is the lattice spacing. Gauge invariance 
implies that the fermion determinant depends only on TrP = P r + iP{. 
It can be calculated by Fourier transform on the lattice, using staggered 
quarks. As this work was in progress, the preprint || appeared, where the 
same calculation of the determinant is performed in the continuum. This 
determinant represents the fermion contribution to the one-loop effective 
potential of the Polyakov loop. It should become a better approximation 
as the gauge fluctuations are suppressed, i.e. at high temperature. The be- 
havior of the determinant as a function of (P r ,Pi) is summarized in Figs. 3 
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to 5 for m q = 0. Increasing the quark mass does not change the qualitative 
behavior, but only reduces the amount of variation of the determinant. 
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Figure 3: Isocontour lines for 1/Vlog | det | in the toy model. 

Fig. 3 shows y log | det |, as a function of the Polyakov loop. It is max- 
imum in the free field case Aq = 0,P = (3,0), in accord with the theorem 
of E. Seiler Q. Note that isocontour lines are almost vertical and almost 
equally spaced, which implies that log | det | is almost a linear function of P r 
alone: y log | det | ~ c(/i, T, m q )P r + cq. What is particularly interesting is 
the variation of the slope c(n, T, m q ) with fi and T, since this quantity char- 
acterizes the amount by which larger Polyakov loop values P r are favored, 
i.e. the ordering effect of the fermion determinant. Fig. 4 shows clearly that 
a similar ordering effect can be achieved by increasing the temperature or 
the chemical potential. 

The oscillatory behavior of the determinant is displayed in Fig. 5, which 
shows the phase y</>(det), as a function of the Polyakov loop. The isocontour 
lines correspond to increases of the phase by 2ir. They are almost parallel 
and equally spaced, confirming <p(det) oc Pi in this model. This figure makes 
it clear how the sign problem occurs: as V increases, the quenched measure 
exp(— S g ) allows for fluctuations in Pj of 0(l/\/V), which are sufficient to 
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Figure 4: Logarithm of the modulus of the determinant for Pj = in the toy 
model: a decrease in temperature (from middle to top curve) can be nicely 
compensated by an increase in the chemical potential (bottom curve). 
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Figure 5: Lines where </>(det)) = 2kn in the (P r , Pj) plane for the toy model, 
on an 8 3 x 4 lattice. The phase is proportional to the lattice volume, so that 
the lines depicted become denser as the volume increases. 
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rotate the phase </>(det) by 0(W), and drive (cos(c/>(det))) to zero. 

Assume now that the fi = QCD theory admits a deconfinement tran- 
sition at T = To, characterized by 2 degenerate minima of the free energy 
at neighboring (or degenerate) real values Pi < P2 of the Polyakov loop. 
If the chemical potential [i is turned on, the fermionic determinant will fa- 
vor the minimum at P2 because of the ordering effect of [i (Fig. 4). The 
balance between Pi and P2 can be restored by a compensating decrease in 
temperature. This is achieved when the derivative -%p-log(\ det(fx, T)\)\p 1 
keeps its (fj, = 0, T = Tq) value. In this way, a heuristic phase boundary 
/i c (T) can be obtained, as in Fig.6. The phase boundary can be obtained 
by the same arguments using the continuum expression for the fermion de- 
terminant in a gauge field background Aq (Eq. [12] of Ref. ||). The lines 
of Fig.6 come out from the continuum expression, and the discrete points 
from the lattice expression for staggered fermions. The two sets of points 
correspond to Pi = 3 and Pi = 0, to give some measure of our systematic 
uncertainties. The phase diagram shows a quadratic dependence of the crit- 
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Figure 6: Heuristic phase boundary obtained from the toy model. The 
lines are derived from the continuum calculation of [||, the data points are 
obtained from the lattice. The 2 boundaries correspond to different critical 
values of P r . 
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Figure 7: Spectrum of the fermion matrix for (P r ,Pj) = (0,0) on a 8 3 x 4 
lattice, for three values of the chemical potential [i. Eigenvalues are repelled 
from the origin as [i increases, indicating chiral symmetry restoration. 



ical temperature versus \x. It even is quantitatively plausible. Nonetheless, 
several crude approximations have been made: all fluctuations of the gauge 
field have been ignored except the zero-mode of the Polyakov loop; and the 
variation of the quenched measure with T is neglected. The failure of this 
crude approximation becomes obvious as the temperature is lowered: the 
critical chemical potential /U C (T) diverges as T — > []. Therefore, it is likely 
that jU c (T) consistently overestimates the true critical \x. 

However, our toy model may still provide some qualitative insight as 
to what happens at low temperatures. As T — > 0, the quenched measure 
for the Polyakov loop becomes Gaussian around the origin: Prob(P) oc 
exp(— cV(Pr + Pi))- On the other hand, note the slight curvature of the 
contour lines in Fig. 3. They indicate that, if P r fluctuates little around 0, 

1 This is brought about by the breakdown of the approximation for the fermion deter- 
minant, by which log(|det(/i = 0)|) varies as T 3 (see ||) instead of 1/T as it should as 
T -» 0. 
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Zo<7(|det|) increases with |Pj|: the determinant favors large values of Pj, 
whereas the quenched measure favors small values. The partition function 
has the form 

JdP ie - cVP ? e wrm+vdMP? (3) 

with the first factor coming from the quenched measure, the second one 
from the fermion determinant. If for some value of fi, d{y) = c, then a 
phase transition ensues. The major contribution to the partition function, 
instead of coming from the region Pj ~ 0, will come from the large values 
of Pi. To support this plausible scenario, Fig. 7 shows the Dirac spectrum 
in the complex plane for our toy model, at (P r ,Pj) = (0,0), for increasing 
values of /x. As \x increases, eigenvalues are repelled from the origin, which 
would be consistent with (if)ip \ = and chiral symmetry restoration at large 
p. 

Finally, a remarkable feature of our heuristic phase diagram is that its 
form remains unchanged for SU(2). There is no difference at all if one 
chooses for Pi its free field value (upper curve in Fig. 6). For any other 
value the difference between the phase boundary for the SU (2) and SU (3) 
cases is a constant factor. An exact SU{2) phase diagram should soon be 
available for comparison, since SU(2) simulations do not suffer from the sign 
problem. 

4 Constrained Monte Carlo 

To go beyond our toy model, in SU(3), one must address the sign problem. 
We have considered the introduction of a constraint on Pj in the partition 
function. Namely we want to study by Monte Carlo the modified partition 
function 

Z'=f [dU] e~ s ^ u) det 2 (f/,/x) 5(Pi(U)) (4) 

The benefits are clear: as seen in Section I, the fluctuations in 0(det) will be 
reduced, especially for heavy quarks. In the same way that the microcanon- 
ical and canonical partition functions become equivalent in the thermody- 
namic limit, the constraint on Pj becomes irrelevant as V —* oo, provided Pj 
is frozen to its mean value given by the saddle point of its effective potential. 
However, at finite density, the mean value of Pj is imaginary, as can be seen 
from Eq.(^|). Therefore, our constraint introduces some systematic error. 
This can be traced to the following. 
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Figure 8: Average of cos(0(det)) as a function of the fermion mass m in 
two quenched simulations on a 4 3 x 6 lattice at (3 = 1 and /i = 1.5. Both 
simulations have equal statistics. One of them includes the constraint Pj = 0. 



The free energy F q of a single quark at finite /i is given by the Polyakov 
loop: 

e -(3F q = / Tr p^ =Pr + jp. ( 5 ) 

where Pj is pure imaginary. On the other hand, for an antiquark one has 

e -pF, = / Xr pt) = Pr _ iP . ( 6 ) 

Constraining Pj = thus enforces equality of the 2 free energies. This is a 
good approximation in 2 regimes: 

- \P r \ <C | P; |, i.e. high temperature or small chemical potential; 

- |P r | = | Pj | =0, i.e. zero temperature, for any /i < /i c . 

This second regime is particularly interesting and difficult to investigate. 

To test our ideas, we have implemented a constrained Hybrid Monte 
Carlo algorithm, and performed quenched 5^7(3) simulations, measuring the 
average "sign" (cos((/>(det))) with and without the constraint Pi = 0. Fig.8 
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shows our results as a function of the quark mass, for chemical potential 
fj, = 1.5 in lattice units, at T ~ 0. The relative error on (cos(</>(det))), which 
propagates to all observables, is reduced by over an order of magnitude as 
one approaches the phase transition. 



5 Conclusion 

We have exhibited the strong correlation present in static QCD between the 
phase of the fermion determinant and the imaginary part of the Polyakov 
loop. This correlation motivated us to study a toy model which keeps the 
zero-mode of the Polyakov loop as the only degree of freedom. This toy 
model clearly shows the ordering effect of the chemical potential and the 
oscillations of the fermionic measure. A heuristic phase diagram can be 
obtained, and a mechanism for a T = phase transition suggests itself. 
Finally, constrained Monte Carlo simulations where the imaginary part of 
the Polyakov loop is fixed to zero provide an important reduction of the sign 
problem, and allow to go beyond static QCD. 

We thank the authors of Ref.[[l|], especially Doug Toussaint, for mak- 
ing their data available to us. Ph. de F. thanks Mike Creutz for helpful 
discussions. 
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